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NOMENCLATURE 


Hertzian  contact  radius 

Modulus  o£  elasticity  of  the  contact  materials  1  and  2 
respectively 

Equivalent  modulus  of  elasticity  of  the  contact  materials, 


where,  E, 


rl 


1 

TT 


Material  parameter,  where  G  =  (aE)  dimensionless 
Elastic  shear  modulus  of  the  lubricant 
film  thickness 
nominal  film  thickness 

film  thickness  parameter,  where  H  =  h/R,  dimensionless 

nominal  film  thickness  parameter,  \diere  ^imensionle 

Length  of  the  contacting  cylinders 

slope  of  the  limiting  shear  stress -pressure  relation 

pressure 

pressure  gradient,  p'  =  , 

volumetric  flow  rate 

radius  of  the  contacting  cylinders 

equivalent  radius  of  the  contacting  cylinders 

where  1/R  =  1/r^  +  l/r2 

time 

lower  surface  velocity  in  direction 
upper  surface  velocity  in  direction 
rolling  speed,  where  u  =  u^^  +  U22 


V 


,  dimensionless 


U 

w 

w/L 

W 


a 

U 


=  speed  parameter  where  U  = 

=  load 

=  load  per  unit  length 
=  load  parameter,  where  W  = 

=  horizontal  coordinate  axis  (origin  is  at  the  center 
o£  the  Hertzian  contact) 

=  vertical  coordinate  axis 
=  pressure  coefficient  of  viscosity 
=  viscosity 
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=  zero  pressure  value  of  viscosity 


^1,2 

^1 

'^2 
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=  Poisson's  ratio  for  contact  materials  1  and  2  respectively 
=  limiting  shear  stress 
=  lower  surface  shear  stress 
=  upper  surface  shear  stress 
=  zero  pressure  value  of  limiting  shear  stress 


=  slide/roll  ratio,  where  Z 

=  — ,  i  =  1  or  2 
'^L 

_  y  dx 

"  TiTT  3t 

L  oo 


2  ^21  '  ^11 
U21  +  u^i 
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I.  INTRODUCTION 


In  a  previous  report  [1]  the  authors  have  shown  that  under  well 
defined  conditions  many  liquid  lubricants  behave  as  non-equilibrium 
amorphous  solids  in  the  Hertzian  region  of  elastohydrodynamic  contacts. 
They  have  also  presented  a  simple  shear  rheological  constitutive 
equation  [2]  which  requires  only  three  material  properties  as  functions 
of  pressure  and  temperature.  These  are  the  low  shear  stress  viscosity 
(y^(p,T)),  the  limiting  elastic  shear  modulus  (G^(p,T)),  and  the  limiting 
shear  stress  (tj^(p,T)). 

The  low  shear  stress  viscosity  of  many  lubricants  has  been 
reported  in  the  literature  for  nearly  a  century  and  is  therefore  a 
familiar  concept.  The  limiting  elastic  shear  modulus  (or  high  frequency 
shear  modulus)  is  less  familiar  in  the  lubrication  literature  but 
nevertheless  has  been  measured  for  many  years  by  several  techniques 
including  ultrasonics  [3].  The  limiting  shear  stress  of  liquid 
lubricants  has  been  the  subject  of  speculation  for  many  years  [4] 
and  indeed  the  nature  of  traction  measurements  in  EHD  contacts 
has  led  several  researchers  to  support  that  view  for  sometime  [5, 6, 7, 8]. 

The  first  measurements  of  limiting  shear  in  liquid  lubricants, 
independent  of  EHD-type  experiments,  were  reported  by  Bair  and  Winer  [9] 
at  temperatures  and  pressures  typical  of  the  Hertzian  region  of  EHD 
contacts.  Because  at  those  temperature -pressure  conditions  the  material 
was  in  the  amorphous  solid  state,  some  if  not  all  of  its  properties 
could  be  expected  to  be  a  function  of  the  history  (state  path)  to  which 
the  material  was  subjected  as  it  went  from  the  liquid  to  the  amorphous 
solid  state  [10].  Qualitatively,  where  the  material  passes  through  the 
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liquid-solid  transition  will  be  one  of  the  defining  characteristics 
of  the  material  in  the  amorphous  state  [10].  The  higher  the  pressure 
at  which  the  transition  occurs,  the  higher  will  be  the  density  and 
limiting  shear  stress  in  the  amorphous  state. 

In  the  previously  reported  limiting  shear  stress  research  [2] 
the  material  had  been  subjected  to  an  isobaric  cooling  process  for 
reasons  of  convenience  in  the  experimental  technique .  However ,  in  a 
typical  EHD  contact  at  low  slide -roll  ratio  the  temperature  rise  is 
small  [11]  and  the  process  the  lubricant  is  subjected  to  as  it  passes 
through  the  contact  is  closer  to  isothermal  compression.  Therefore, 
for  the  same  temperature  and  pressure,  the  material  in  the  limiting 
shear  stress  experiment  had  passed  through  the  liquid-solid  transition 
at  a  higher  pressure  (therefore  higher  density  in  the  amorphous  solid) 
than  would  occur  in  the  EHD  contact.  This  might  be  expected  to  cause  the 
measured  limiting  shear  stress  to  be  higher  (at  least  different)  than 
would  be  the  case  in  the  EHD  contact  at  the  same  temperature  and  pressure. 
This  may  also  explain  why  the  predicted  traction-slide-roll  ratio  curves 
[2]  based  on  the  laboratory  limiting  shear  stress  measurements  had  a 
maximum  about  ten  percent  higher  than  the  EHD  measured  values  even 
though  the  increasing  portion  of  the  curves  agree  quite  well .  One  of 
the  objectives  of  this  research  was  to  examine  the  effect  of  history 
on  the  limiting  shear  stress  of  liquid  lubricants.  It  will  be  shown 
that  the  isothermal  compression  history  results  in  a  lower  limiting 
shear  stress  and  one  which  agrees  with  the  peak  traction  measurement 


in  EHD  contacts. 
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The  shear  stress -strain  behavior,  including  limiting  shear 
stress,  of  solid  lubricating  plastic  materials  have  also  been  measured 
under  pressure.  The  history  of  the  solid  plastic  lubricating  materials 
in  a  lubricated  contact  would  not  be  expected  to  be  important  in  deter¬ 
mining  the  shear  properties  of  the  lubricant  because  they  do  not  experi¬ 
ence  a  liquid-solid  transition  in  the  contact.  However,  the  previous 
processing  history  of  the  material,  when  it  may  have  passed  between 
liquid  and  solid  states,  will  influence  the  properties.  The  previous 
processing  history  is  known  to  influence  the  yield  stress  of  solid 
polymers  [12]. 

The  shear  rheological  model  with  measured  limiting  shear  stress 
was  also  used  to  predict  EHD  traction  slip  roll  ratio  behavior.  Those 
predictions  were  compared  with  EHD  data  published  by  Johnson  and 
Tevaarwerk  and  measured  in  our  laboratory.  Side  slip  and  spin  were 
introduced  into  the  model.  Side  slip  due  to  misalignment  was  shown  to 
be  very  influential  on  traction  both  in  the  model  and  the  measurements. 

A  Grubin-like  EHD  inlet  analysis  utilizing  a  non-linear  viscous 
fluid  model  with  a  limiting  shear  stress  is  also  reported.  The  shear 
rheological  equation  requires  only  a  low  shear  stress  viscosity  and 
the  limiting  shear  stress  both  functions  of  pressure.  Values  employed 
for  these  properties  are  taken  from  measurements  on  typical  lubricants . 
Reductions  of  EHD  film  thickness  are  found  to  be  up  to  40  percent 
compared  with  the  standard  Grubin  prediction  for  typical  operating 
conditions.  Slide-roll  ratio,  limiting  shear  stress  dependence  on 
pressure,  and  atmospheric  pressure  value  of  limiting  shear  stress  are 


new  variables  required  to  determine  film  thickness  with  the  first  two 
being  more  important  than  the  last.  The  EHD  film  thickness  is  reduced 
by  increasing  slide-roll  ratio  and/or  decreasing  the  pressure  dependence 
of  the  limiting  shear  stress. 
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II.  EXPERIMENTAL  APPARATUS 

Two  apparatus  were  used  to  measure  the  high-pressure  shear 
rheological  response  o£  lubricants  to  1.1  GPa  (162,000  psi) .  One, 
the  low  pressure  stress- strain  apparatus  was  used  in  some  of  the 
isobarically  cooled  investigations  and  the  other,  referred  to  as  the 
constant  pressure  stress- strain  apparatus,  was  used  in  all  of  the 
isothermal  compression  experiments  and  all  of  the  experiments  involving 
solid  polymers. 

A.  Low  Pressure  Stress -Strain  Apparatus:  0.7  GPa 

As  reported  previously  [2],  an  apparatus  was  constructed  to 

measure  the  mechanical  shear  properties  of  glassy  lubricant  samples 
to  pressures  of  0.7  GPa.  It  is  shown  schematically  in  Figure  1.  The 
glassy  sample  is  formed  in  an  annular  groove  by  cooling  at  elevated 
pressure.  The  groove  is  kept  filled  by  a  sample  reservior  which  is 
sealed  from  the  working  fluid  (gasoline)  by  an  isolator  piston.  The 
sample  material  can  be  sheared  in  the  annulus  by  the  development  of 
a  pressure  difference  across  the  driving  piston.  The  shear  stress  is 
determined  by  knowing  the  geometry  and  measuring  the  differential 
pressure  by  two  pressure  transducers.  The  sample  strain  is  determined 
by  the  displacement  of  the  driving  piston  measured  with  an  LVDT.  This 
signal  can  also  be  used  to  measure  the  strain  rate.  The  shear  stress 
(pressure  difference) ,  the  strain  (piston  displacement)  and  time  are 
recorded  on  an  x-y-y  recorder.  Sample  temperature  is  determined  by  a 
thermocouple  imbedded  in  the  pressure  vessel  wall. 

At  moderate  working  temperatures  such  as  those  for  5P4E  (-20  to  35C) , 
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and  elevated  pressures,  the  seal  friction  is  negligible  and  no  shearing 
force  across  the  piston  can  be  maintained  when  the  test  material  is 
above  its  glass  transition  temperature.  However,  with  Nl*  the  tempera¬ 
ture  required  to  go  into  the  glassy  region  at  moderate  pressures  is  so 
low  (-40C)  that  a  correction  for  seal  friction  must  be  employed.  The 
seal  friction  at  low  temperature  was  calibrated  by  using  gasoline  as 
the  test  fluid  which  has  very  low  viscosity  at  the  test  temperature  and 
pressure.  Therefore,  at  the  low  shearing  rate  of  the  experiment,  the 
driving  force  on  the  piston  was  assumed  to  be  due  to  seal  friction.  This 
seal  friction  was  typically  less  than  five  percent  of  the  maximum  shear 
stress  measured  for  N1 . 

Referring  to  Figure  1  the  sequence  of  a  typical  experiment  is 
the  following;  with  the  sample  in  the  apparatus,  the  system  is  heated 
to  a  temperature  high  enough  to  keep  the  sample  in  the  liquid  region 
at  the  predetermined  pressure  to  be  used.  The  system  is  then  brought 
up  to  pressure  with  the  valve  open  insuring  uniform  pressure  through¬ 
out  the  apparatus.  The  system  is  then  cooled  to  the  desired  tempera¬ 
ture  at  or  below  the  dilatometric  liquid-solid  transition  while  main¬ 
taining  constant  pressure.  The  isolating  piston  movement  accommodates 
sample  volume  change  during  these  state  changes .  The  valve  is  then 
closed  isolating  the  regions  above  and  below  the  driving  piston.  Stress 
is  applied  to  the  sample  by  either  increasing  or  decreasing  the  pressure 
on  the  bottom  of  the  driving  piston  by  varying  the  supply  pressure. 

The  pressure  difference  is  measured  by  the  two  pressure  transducers . 


*The  fluids  are  described  in  the  Appendix. 
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'rhe  driving  piston  displacement  and  velocity  are  measured  by  the  LVDT. 

By  the  nature  o£  the  device,  when  the  piston  moves  downward  the  pressure 
level  decreases  and  when  it  moves  upward  the  pressure  level  increases . 
The  pressure  level  changes  can  be  kept  to  a  minimum  by  keeping  the 
strain  (piston  displacement)  small  for  a  given  measurement. 

B.  Constant  Pressure  Stress -Strain  Apparatus  1.1  GPa 

Previous  shear  strength  measurements  [1,2]  were  made  in  this 
laboratory  with  an  apparatus  which  derived  the  force  necessary  to 
shear  the  sample  from  the  pressure  of  the  pressurizing  medium.  Conse¬ 
quently,  as  the  sample  was  strained  the  hydrostatic  pressure  changed. 

A  new  apparatus.  Figure  2,  has  been  constructed  to  perform  at  a  nearly 
constant  pressure,  extend  the  pressure  range  to  1.1  GPa  (162  kpsi) ,  and 
accommodate  large  changes  in  volume  of  liquid  sarples.  In  addition,  a 
replacement  cell  for  the  testing  of  solid  polymer  samples  has  been 
provided. 

The  apparatus  includes  an  integral  pressure  intensifier  whose  high 
pressure  piston  forms  one  of  the  closures  of  a  translating  cylindrical 
pressure  vessel.  The  other  closure  is  a  fixed  piston.  The  vessel  can 
be  driven  hydraulically  by  oil  supplied  to  either  end  of  the  vessel. 

For  solid  polymers  the  experimental  cell  enclosed  in  the  vessel  is 
simply  a  fixture  for  holding  and  shearing  an  annulus  shaped  sample. 

For  liquid  samples,  the  cell  includes  a  reservoir  with  an  isolating 
piston  to  replenish  the  sample  in  the  annulus  as  its  volume  is  reduced 
on  pressurization. 

The  vessel  can  be  moved  a  small  distance  before  the  cell 
contacts  the  fixed  piston  allowing  calibration  of  the  closure  seal 
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friction.  Further  movement  of  the  vessel  produces  a  relative  motion 
between  the  cell  and  the  vessel,  straining  the  sample.  The  vessel 
displacement,  hence  sample  strain,  is  measured  by  an  LVDT  transducer. 
The  hydraulic  pressure  which  drives  the  vessel  is  measured  by  a 
commercial  pressure  transducer.  The  strain  rate  is  controlled  by 
regulating  the  hydraulic  oil  flow  rate.  Temperature  is  determined 
with  a  copper-constantan  thermocouple  in  the  wall  of  the  vessel. 

The  apparatus  is  submerged  in  an  oil  bath  for  temperature  control. 

The  pressure  of  the  large  end  of  the  intensifier  is  measured 
with  a  precision  Heise  bourdon  tube  gauge  and  the  pressure  of  the 
medium  in  the  vessel  (diester)  is  determined  from  the  known  intensifier 
area  ratio  and  measured  seal  friction.  Hydraulic  pressure  acting  on 
the  vessel  (sample  stress)  and  vessel  displacement  (sample  strain)  are 
recorded  as  functions  of  time  on  a  two  function  x-y-y  plotter. 


III.  EXPERIMENTAL  LUBRICANTS 


The  materials  included  in  this  study  were 


Liquids 

5P4E 

N1 

N1  +  Polymer 
MCS  460 
KRYTOX 
Santotrac  50 
XRM  177F 
LVI  260 
Vitrea  79* 


Solids 

Teflon 

Polyvinylchloride 
Acrylic  (extruded) 
Nylon 


Description  of  liquids  may  be  found  in  Appendix  A. 
polymers  were  obtained  form  a  local  plastic  supply  house 
and  machined  to  fit  the  apparatus. 


Solid 

in  rod  form 


*Sanples  received  from  K.  L.  Johnson,  Cambridge  University. 


10 


IV.  EXPERIMENTAL  RESULTS 

The  data  take  the  form  of  stress  versus  strain  curves  obtained 
from  an  x-y  plotter  along  with  a  plot  of  time  versus  strain  from 
which  strain  rate  is  determined.  Experiments  involving  isobaric  cooling 
history  were  begun  at  atmospheric  pressure  above  the  glass  transition 
temperature  as  determined  by  dilatometry,  pressurized  at  that  temperature, 
then  cooled  to  the  test  temperature.  To  determine  limiting  shear  stress 
(t,  )  each  sample  was  sheared  plastically  at  various  strain  rates. 

JL 

Figure  3  is  a  typical  plot  of  shear  stress  versus  shear  strain-rate 

for  the  lubricant  Nl.  At  the  lowest  pressure,  viscous  behavior  is 

exhibited  at  the  lower  strain  rates  by  the  slope  of  the  curve  tending 
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toward  45  degrees.  The  viscosity  of  the  sample  must  be  at  least  10  Pas 
at  temperature  and  pressure  in  order  to  measure  limiting  shear  stress 
due  to  the  limited  shear  strain  rate  provided  by  this  instrument. 

A.  Liquid  Lubricants 
1.  Limiting  Shear  Stress 

As  a  continuation  of  previous  work  [2],  limiting  shear  stress  was 
measured  for  Krytox  (perfluorinated  polyether)  and  Nl  plus  four  percent 
PAMA  (polyalkylmetacrylate)  at  one  pressure  and  MCS  460  (synthetic 
hydrocarbon)  at  two  pressures ,  (Figure  4) .  Nl  (naphthenic  base  oil) 
is  also  shoivn  for  comparison.  The  history  is  isobaric  cooling.  A 
small  reduction  in  Tj^  results  from  the  addition  of  polymer  to  Nl , 
consequently  some  loss  in  EHD  traction  would  be  expected. 

Since  it  is  known  [10]  that  the  history  of  an  amorphous  material 
as  it  goes  from  the  liquid  to  the  solid  state  influences  the  resulting 
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density  o£  the  amorphous  solid,  an  experiment  was  performed  with  5P4E 
to  determine  the  effect  of  history  on  limiting  shear  stress.  The  upper 
curve  in  Figure  5  represents  the  limiting  shear  stress  versus  pressure 
at  38C  resiilting  from  isobaric  cooling  in  the  constant  pressure  shear 
stress  apparatus .  Included  are  data  points  from  the  low  pressure 
apparatus  and  high  stress  viscometer  previously  reported  [2].  The 
curve  below  represents  data  for  the  same  material  and  temperature 
subjected  to  isothermal  compression.  The  compression  history  produces 
a  lower  than  the  cooling  history  as  would  be  expected  since  the 
density  is  also  lower  for  the  compressional  history.  Also  in  Figure  5 
is  the  limiting  shear  stress  for  5P4E  at  80C  for  compressional  history. 
Since  isothermal  conpression  is  the  history  most  representative  of  the 
inlet  of  an  EHD  zone  where  a  solid  transition  may  occur,  all  succeeding 
data  will  be  for  that  history. 

In  Figure  6,  pres sure -limiting  shear  stress  isotherms  are  plotted 

for  N1  (naphthenic  base  oil) ,  Krytox  (perfluorinated  polyether) , 

Santotrac  50  (cycloaliphatic  hydrocarbon  traction  fluid) ,  and  the  samples 

of  K.  L.  Johnson,  Vitrea  79  and  LVI  260.  An  attempt  was  made  to  measure 

the  stress -strain  behavior  of  XRM  177F  to  1.1  GPa  at  22C,  but  no 

measurable  shear  stress  was  developed  at  the  available  shear  rate  because 

2 

of  its  low  viscosity  (<  10  Pas)  at  that  pressure  and  temperature.  It 
is  noteworthy  that  the  order  of  the  magnitude  of  for  Santotrac  50,  N1 , 
and  Vitrea  79  is  the  same  order  in  which  one  would  expect  to  rank  them 
for  traction  coefficient. 
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2.  Shear  Modulus 

The  slope  o£  the  initial  linear  portion  o£  the  shear  stress - 
shear  strain  curve  yields  the  elastic  shear  modulus  o£  the  lubricant 
sample  when  corrected  for  the  deflection  of  the  instrument  parts  during 
the  experiment.  The  instrument  deflection  was  calibrated  by  replacing 
the  sample  cell  by  the  solid  steel  slug.  Figure  7  shows  the  shear 
modulus,  G,  of  N1  versus  shear  rate  at  four  pressures.  At  0.91  GPa 
pressure  the  shear  modulus  is  independent  of  strain  rate  for  the  range 
of  rates  shown.  It  is  expected  that  at  the  lower  pressures  a  limiting 
shear  modulus,  G^,  would  also  be  reached  if  the  strain  rate  were  high 
enough  or  the  temperature  low  enough. 

B.  Solid  Polymers 

The  shear  rheology  of  polyvinyl  chloride.  Teflon,  acrylic,  and 
nylon  are  shown  in  Figure  8-11.  Due  to  the  volume  contraction  of  the 
sanple  under  pressure  some  clearance  occurs  between  the  polymer  sample 
and  the  holding  fixture  which  prevents  an  accurate  determination  of 
shear  strain.  Therefore,  the  abscissa  of  Figures  8-11  is  relative 
deflection  of  the  cell  components.  This  deflection  is  proportional  to 
sample  strain  in  such  a  ratio  that  1  mm  deflection  is  a  strain  of 
approximately  120  percent.  The  area  used  in  the  shear  stress  calcula¬ 
tion  is  the  area  of  the  sample  at  atmospheric  pressure.  The  rate  of 
deflection  for  the  solid  polymer  measurements  was  approximately  constant 
at  0.2  mm/s  which  is  a  strain  rate  of  approximately  0.24  s 

Since  the  viscosity  of  the  polymer  samples  must  be  very  high, 
the  visco-plastic  relaxation  time  [2] 
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must  be  long  enough  to  produce  the  limiting  shear  stress  at  the  rate 
encountered  in  these  experiments.  The  ultimate  or  limiting  shear 
stress  for  PVC  and  Teflon  is  plotted  in  Figure  12  along  with  the 
limiting  stress  of  N1  and  Santotrac  50.  Not  only  are  the  magnitudes 
comparable  but  the  slopes,  are  approximately  the  same. 
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V.  APPLICATION  OF  THE  MODEL  TO  HID  TRACTION  PREDICTIONS 

The  above  mentioned  isothermal  compression  limiting  shear  stress 
and  limiting  elastic  shear  modulus  were  utilized  in  the  constitutive 
equation  developed  in  the  previous  report  [2]  along  with  the  low  rate 
viscosity  ]j^(P,T]  to  predict  EHD  traction.  The  contact  was  divided 
nonuniformly  into  a  grid  o£  20  segments  on  a  cord  in  the  direction  of 
motion  and  20  such  strips  across  the  contact  perpendicular  to  the 
direction  of  motion  to  permit  pressure  and  material  property  variation 
in  the  contact.  The  following  assumptions  were  employed;  the  film 
thickness  and  material  temperatures  were  assumed  uniform  throughout 
the  contact,  the  pressure  distribution  was  Hertzian,  the  viscosity 
was  an  exponential  function  of  pressure,  the  elastic  shear  modulus 
was  proportional  to  pressure,  the  elastic  surface  compliance  was 
proportional  to  the  contact  traction  as  developed  by  Kalker  [13]  and 
reported  in  Johnson  and  Roberts  [6],  and  inlet  zone  effects  were 
neglected.  Several  of  these  assumptions  can  be  called  into  question 
and  should  be  refined  in  subsequent  development  particularly  those 
concerned  with  the  temperature  distribution  and  the  inlet  zone  influence. 

Although  we  know  the  film  temperature  is  not  constant  the 
analysis  is  done  for  slide-roll  ratios  of  less  than  one  tenth.  From 
other  work  in  this  laborator>^  under  conditions  similar  to  those  used 
in  this  analysis  we  know  the  maximum  surface  temperature  rise  is  usually 
less  than  5C  above  the  bulk  temperature  in  this  range  of  operating 
conditions.  Although  we  have  not  measured  lubricant  temperatures  at 
these  low  slide-roll  ratios,  work  in  sliding  contacts  would  indicate 
they  are  probably  less  than  5  to  IOC  above  the  surface  temperature. 
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With  the  above  assumptions  a  program  was  written  to  calculate 
the  local  shear  stress  at  each  point  in  the  grid  by  using  a  Bisection 
Method  on  the  model  equation  with  starting  shear  stresses  o£  zero 
and  0.999  x.  .  I£  the  Bisection  Method  does  not  £ind  a  solution  as  the 

Li 

trial  shear  stress  reaches  0.999  x^^,  the  solution  is  assumed  to  be 
To  obtain  the  time  derivative  term,  upstream  grid  positions  plus  a 
convective  derivative  are  employed  £or  a  given  grid  point.  The  average 
shear  stress  in  the  contact  is  obtained  by  integration  over  the  area 
and  the  traction  coe££icient  is  the  ratio  o£  the  average  shear  stress 
divided  by  the  average  pressure. 

Figure  13  shows  the  three  predicted  traction  curves  £or  5P4E, 

LVI  260,  and  Vitrea  79  at  an  average  Hertzian  pressure  o£  0.67  GPa 
and  the  indicated  temperatures  assuming  that  no  spin  or  side-slip  was 
present.  For  comparison  the  experimental  data  reported  by  Johnson  and 
Tevaarwerk  [8]  are  included.  The  prediction  o£  peak  traction  seems  to 
be  in  good  agreement  with  experiment,  however,  the  low  slide -roll  ratio 
portiono£  the  predicted  curves  lie  at  about  one  third  the  slide -roll 
ratio  o£  the  experimental  points. 

In  order  to  reconcile  the  di££erence  at  low  slide-roll  ratio  o£ 
the  predicted  curves  and  experiment  both  side-slip  and  spin  were  added 
to  the  traction  program.  They  can  be  added  easily  as  two-dimensional 
strain  to  the  visco-plastic  or  elastic  plastic  equations  [2],  but  a 
£ull  visco-elastic-plastic  solution  £or  spin  and  slip  has  not  been 
completed.  I£,  however,  the  viscosity  is  very  large,  then  the  viscous 
portion  o£  the  solution  may  be  omitted  by  setting  very  large  without 
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altering  the  predicted  traction  curve.  In  the  same  manner,  if  the 
viscosity  is  small,  then  elasticity  may  be  omitted.  The  spin  and  side 
slip  were  considered  because  they  may  have  been  present  in  the  experi¬ 
ment  as  a  result  of  small  misalignment  and  will  be  present  to  some 
extent  in  most  EHD  applications . 

In  Figure  14  side  slip  has  been  added  to  the  predictions  for 
5P4E  at  two  pressures.  When  comparing  these  predictions  with  the 
Johnson  and  Tevaarwerk  data  [8]  a  side  slip  angle  (as  defined  by  the 
illustration)  of  only  0.6  milli-radian  was  required  to  bring  the 
low  slide-roll  ratio  portion  of  the  model  into  agreement  with  experi¬ 
ment  as  shown  in  Figure  15.  It  was  assumed  that  no  spin  was  present. 

A  small  amount  of  spin  or  combination  of  slip  and  spin  would  have  the 
same  effect.  Such  a  small  angular  misalignment  of  machine  elements 
is  not  inprobable  in  even  the  most  carefully  assembled  EHD  simulators. 

To  further  evaluate  the  effect  of  small  side  slip  resulting  from 
axis  misalignment  our  ball  on  flat  HID  simulator  was  modified.  The 
modifications  considered  of  an  ability  to  adjust  the  sapphire  mount 
to  permit  movement  of  the  sapphire  axis  of  rotation  relative  to  the  ball 
axis  of  rotation.  Although  precise  alignment  could  not  be  determined, 
the  relative  angle  could  be  varied  by  hundreds  of  a  radian.  The  angle 
of  interest  is  that  formed  between  two  lines  in  the  plane  of  the 
sapphire  surface  both  passing  through  the  hertz  contact,  one  parallel 
to  the  ball  axis  of  rotation  and  one  through  the  point  where  the  sapphire 
axis  of  rotation  passes  through  the  surface.  Precise  alignment  occurs 
when  the  sapphire  axis  intersects  the  ball  axis  of  rotation.  An 
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additional  modification  was  the  ability  to  continuously  record  surface 
velocities  and  traction.  These  signals  were  digitized  and  analyzed  with 
a  program  to  determine  the  traction-slide/roll  ratio  curves  while  the 
surface  velocities  were  driven  in  such  a  manner  as  to  sweep  the  slide 
roll  ratio  from  negative  to  positive. 

Figure  16  shows  the  traction  coefficient  versus  slide -roll  ratio 
relation  for  the  mineral  oil  N1  at  the  conditions  indicated  and  several 
side  slip  angles  as  predicted  by  the  model,  Figure  17  contains  traction 
curves  measured  for  N1  and  the  same  conditions  and  four  relative  side 
slip  angles.  By  comparing  the  predicted  (Figure  16)  and  measured 
(Figure  17)  we  conclude  the  zero  relative  side  slip  in  the  measured  data 
could  be  0.015  radians.  The  predicted  values  from  the  model  for  the 
same  measured  relative  side  slip  are  shown  in  Figure  18.  The  agreement 
between  measured  and  predicted  values  is  quite  good  ’uhich  lends 
credability  to  the  rheological  model  and  the  limiting  shear  stress 
concept . 

Figure  19  shows  a  measured  traction  curve  for  Santotrac  50.  The 
side  slip  and  spin  were  minimized  within  the  present  limitations  of 
the  apparatus .  The  full  model  prediction  was  not  calculated  but  the 
maximum  traction  coefficient  is  in  agreement  with  the  limiting  shear 
stress  property  measurement  for  this  lubricant. 

These  traction  measurement  comparisons  with  the  model  point 
to  the  validity  of  the  limiting  shear  stress  model  and  the  concept  of 
the  limiting  shear  stress  as  an  important  material  property  for  deter¬ 
mining  traction  in  thin  film  lubrication. 
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VI.  APPLICATION  OF  THE  MODEL  TO  E® 

FIIM  THICKNESS  PREDICTION 

The  shear  rheological  model  with  a  limiting  shear  stress  used 
in  the  previous  section  to  determine  contact  traction  has  also  been 
used  in  a  Grubin-type  EHD  film  thickness  analysis.  This  analysis  was 
conducted  by  Mr.  B.  Gecim  as  his  M.S.  (M.E.)  thesis  research  [14]. 

The  study  of  the  mechanics  of  elastohydrodynamic  (EHD)  contacts 
is  primarily  concerned  with  the  film  thickness  developed  and  the  traction 
force  between  the  two  surfaces.  It  is  well  accepted  that  these  two 
aspects  of  the  contact  are,  in  a  sense,  separate  phenomena  in  that  the 
film  thickness  is  determined  by  flow  in  the  inlet  region  and  the  traction 
is  determined  by  phenomena  in  the  Hertzian  region  of  the  contact.  The 
film  thickness  generation,  and  lubricant  behavior  relevant  to  it,  has 
been  better  understood  than  the  traction  behavior. 

In  our  previous  report  [2]  we  have  proposed  a  visco -elastic-plastic 
flow  shear  rheological  constitutive  equation  for  the  lubricant  based  on 
primary  laboratory  property  measurements  which  appears  to  predict  traction 
behavior  in  EHD  contacts.  In  that  model  the  lubricant  reaches  a  limiting 
shear  stress  value  as  a  result  of  a  visco -plastic  or  elastic -plastic 
transition  that  is  related  to  the  glass  transition  under  pressure.  We 
have  also  shown  this  kind  of  behavior  to  occur  in  several  lubricants. 

The  possibility  is  raised  of  the  limiting  shear  stress  behavior  influ¬ 
encing  EHD  film  thickness  generation  if  the  transition  were  to  occur  in 
the  inlet  zone  of  the  contact.  The  objective  of  this  study  was  to  examine 
that  possibility. 

The  visco -plastic  portion  of  the  constitutive  equation  was  modified 
and  coupled  with  the  equations  of  motion  and  conservation  of  mass  and 
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used  in  a  Grubin-type  EHD  inlet  analysis.  Isothermal,  steady, 
incompressible  conditions  with  negligible  body  and  inertial  forces  are 
considered.  The  Grubin-type  film  thickness  analysis  is  performed  for 
an  EHD  line  contact  configuration  in  the  fully-flooded  full -film  regime. 
The  objective  of  this  study  was  to  determine  the  effect  of  lubricant 
limiting  shear  stress  on  the  nominal  film  thickness . 

Two  new  parameters  are  introduced  into  a  classical  Grubin-type 
film  thickness  analysis;  sliding  speed  and  lubricant  limiting  shear 
stress.  The  limiting  shear  stress  is  expressed  as  a  linear  function 
of  pressure,  at  a  constant  temperature.  In  an  EHD  contact  inlet,  pressure 
increases  by  orders  of  magnitude,  hence  the  limiting  shear  stress  may 
also  increase  greatly.  Therefore,  only  under  the  severe  operating 
conditions  of  high  sliding  speeds  or  pressure  gradients ,  does  shear 
stress  reach  the  limiting  value.  Under  such  circumstances  a  decrease 
of  up  to  forty  percent  from  Grubin's  prediction  of  nominal  film 
thickness  was  found. 

Although  the  confirmation  of  the  existence  of  a  limiting  shear 
stress  is  relatively  new  in  the  field  of  tribology,  it  is  an  intrinsic 
material  property  of  the  lubricants.  Therefore  it  is  reasonable  to 
conclude  that,  under  usual  operating  conditions,  with  most  conventional 
lubricants,  no  drastic  changes  or  sudden  collapse  in  film  thickness 
is  expected.  Hence  the  small  ( <  40  percent)  reduction  in  film  thickness 
predicted  in  this  analysis  is  reasonable. 

In  this  analysis  slip  between  the  lubricant  and  the  boundary  was 
not  permitted.  The  usual  no -slip  boundary  conditions  of  viscous  fluid 
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mechanics  was  employed.  An  alternate  approach  to  this  problem  has 
been  suggested  by  Wilson  to  employ  the  solution  technique  used  by 
Aggarwal  and  Wilson  [l^  in  metal  die  drawing  operations  with  a 
lubricant  model  exhibiting  a  limiting  shear  stress.  That  method  was 
also  employed  and  is  reported  as  part  of  reference  P-4] .  As  discussed 
in  rl4|  we  believe  the  Aggarwal  and  Wilson  P-^  model  fails  to  satisfy 
the  physical  constraint  of  conservation  of  mass  when  the  limiting  shear 
stress  is  reached  at  either  surface. 

Bell  [1^,  proposed  an  analytical  solution  model  for  film  thickness 
between  contacting  cylinders,  with  the  Ree -Eying  constitutive  equation 
which  also  has  a  shear  thinning  behavior  but  no  limiting  shear  stress . 

The  analysis  aimed  to  explain  the  phenomena  in  the  case  of  high  rolling 
speeds,  and  high  viscosity  where  Grubin's  prediction  fail  to  match  the 
experimental  observations.  However,  the  analysis,  being  restricted 
to  the  pure  rolling  case,  lost  the  generality  in  comparison  with  the  real 
contact  phenomena,  whereas  the  present  study  analyzes  the  sliding  as  one 
of  the  major  causes  of  film  thickness  reduction.  Bell's  analysis  suffered 
from  the  lack  of  data  on  the  Ree-Eyring  parameters  whereas  the  present 
study  does  have  the  experimentally  determined  lubricant  parameters  (i.e., 
limiting  shear  stress  parameters)  to  put  in  the  analysis. 

In  a  similar  analysis  to  Bell's  earlier  work.  Bell  and  Kannel  P'^ 
analyzed  a  generalized  pressure  dependency  of  viscosity,  non -Newtonian 
Rheology  of  Ree-Eyring  form  and  a  time  delay  in  pressure  effect  on 
viscosity.  The  time  delay  approach  was  found  to  provide  the  best 
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correlation  with  experimental  measurements.  However,  this  analysis 
required  to  specify  an  inlet  pressure  which  is  not  a  well  defined 
parameter  in  the  field. 
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A.  Shear  Rheological  Model 

The  lubricant  shear  rheological  model  employed  is  a  slightly 
modified  version  of  the  limiting  shear  stress  model  put  forth  by  Bair 
and  Winer  [2]  based  on  primary  measurements.  Their  equation  was 


9  u  _  ■  _  J_  dx 
'dx  ^  G  Ht 
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which  can  be  viewed  as  a  modified  Maxwell  viscoelastic  model  with  a 
non-linear  viscous  term  exhibiting  a  limiting  shear  stress,  The 

viscous  term  in  that  model  suffers  from  the  obvious  lack  of  symmetry 
about  zero  shear  stress  and  is  unsuitable  for  incorporating  in  an 
analysis.  To  correct  this  difficulty  it  is  proposed  to  replace  the 
natural  logarithm  with  the  inverse  hyperbolic  tangent  which  has  essen¬ 
tially  the  same  behavior  in  the  positive  shear  stress  range  and  has 
the  advantage  of  having  the  necessary  symmetry.  Therefore,  the 
modified  Maxwell  model  becomes  in  dimensionless  form 


Y  =  T  +  tanh  (t) 

or  in  dimensional  form  Equation  (1)  becomes , 
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From  Equation  (2)  it  is  seen  that  the  three  primary  physical 
properties  required  to  use  the  model  are  low  shear  stress  viscosity 
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[j  the  limiting  elastic  shear  modulus  G^,  and  the  limiting  shear  stress 
T  ,  all  as  functions  of  pressure  and  temperature. 

1j 

For  very  large  values  of  the  limiting  shear  stress  (small  values 
of  fL)  this  model  reduces  to  the  Maxwell  model 


where  the  viscous  term  is  classical  Newtonian  behavior.  The  limiting 
case  of  a  visco -plastic  liquid  results  from  specifying  a  small  visco¬ 
elastic  relaxation  time  ^  ,  so  the  model  takes  on  the  non-linear 

oo 

viscous  form  of 


(2b) 


where 


and 


y  =  yQexp(a-p) 


It  is  this  form  which  is  used  in  this  analysis.  The  magnitude  of  the 
neglected  elastic  term  in  the  resulting  flow  is  small  as  will  be  seen. 


B .  Derivation  of  the  Governing  Equations 
Three  basic  equations  are  used  in  an  EHD  Grubin-type  inlet 
analysis.  They  are  shear  constitutive  equation  (2b),  the  equation  of 
motion  and  conservation  of  mass.  Because  we  assume  isothermal  condi¬ 
tions  throughout,  the  energy  equation  is  not  required.  The  equations 
of  elasticity  are  not  required  because  the  inlet  film  shape  is  assumed 
to  be  known  from  the  Hertzian  contact  analysis . 

The  details  of  the  derivation  are  presented  in  Appendix  B, 
the  resulting  governing  equations  are: 

/T  ^  \ 

^21  ■  ^11  "  ^2^  ■  ^1^ 

+  (1  -  12)  ^nCl  -  T2)  -  (1  -  T^)  £nCl  -  fj^)j  (3-a) 

^  1^(1  +  f2^^^to(l  +  ^2^  -  1/2) 

-  (1  +  ip^(^n(l  +  -  1/2)] 

-  ^  [(1  -  f2)^(5^n(l  -  T2)  -  1/2) 

-  (1  (£n(l  -f^)  -  1/2)]  -  (1  +  x^)h£n(l  +  f^) 


-  (1  -  f3^)hJln(l  -  fp 


(3-b) 
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From  Equation  ^^'^2  ~  ^  F(x^)  -  Tj^(x^).  Since  many  terms 

in  Equations  (3-a)  and  (3-b)  can  be  expressed  as  a  function  of  pressure, 
pressure  gradient  and  x^,  these  two  equations  supposedly  could  be  put  in 
the  form  of  a  simple  first  order  ordinary  differential  equation.  But 
since  it  is  apparently  impossible  to  achieve  this  form  by  making 
elementary  algebraic  manipulations;  equations  (3-a)  and  (3-b)  are  solved 
numerically. 

The  film  thickness  equation  used  in  the  inlet  region  is  h  =  h^  +  h^ 

where 


equation  of  elastic  deformation  outside  the  Hertzian  contact,  due  to 
the  pressure  only  in  the  Hertzian  contact,  which  is,  for  heavily  loaded 
contacts ,  near  Hertz  pressure  distribution  [1^ . 
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•  Method  of  Solution 

The  two  governing  equations  derived  above,  Equations  (3-a) 
and  (3-b),  are  solved  numerically,  for  the  two  unknowns  p'  and 
Since  some  terms  are  functions  of  pressure,  the  solution  method  can 
be  described  as  the  solution  of  two  nonlinear  equations  along  with 
the  solution  of  the  differential  equation  (Equation  5) ,  to  find  the 
pressure  distribution  in  the  inlet  zone. 

The  solution  starts  in  the  inlet  at  where  the  film  thickness 
h  is  ten  times  greater  than  the  nominal  film  thickness  h^ .  Based  on 
starvation  analysis  of  contacts  it  is  reasonable  to  exclude  the 

region  where  h  >  10  h^,  and  perform  the  film  thickness  analysis  in 
the  region  h  <  10  h^  .  The  value  of  x^  at  this  point  is  found  from 
the  film  thickness  equation  (4)  for  each  given  set  of  operating  condi¬ 
tions  and  h^.  At  this  point  pressure  is  assumed  to  be  zero,  and  the 
two  non-linear  equations  are  solved  numerically  for  p'  and  With 

these  initial  conditions  of  p'  and  p.  Equation  (5)  is  solved 
numerically,  to  predict  the  pressure  at  the  next  grid  point.  This 
process  continues  up  to  the  point  where  x^  =  a,  the  Hertz  contact 
radius.  The  distance  between  the  initial  x^  value  and  the  final  value 
of  Xj^  is  uniformly  divided  into  two  hundred  steps  (Ax^) .  Increasing 
the  number  of  grid  points  gives  more  accuracy  but  at  the  same  time 
consumes  more  confuting  time.  Based  on  the  observation  of  a  few  runs 
of  the  program  it  appears  that  the  number  of  grid  points  has  little 
influence  when  increased  from  200  to  500. 
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The  IMSL  library  subroutine  ”ZSYS™"  was  used  for  the  solution 
of  the  two  non-linear  equations.  At  each  grid  point  the  solution  of  the 
previous  point  is  used  as  the  initial  guess  value.  Euler’s  method  is 
used  for  the  solution  of  the  differential  equation 

=  p^  +Ax^  p[  ,  1  =  l,n  (5) 

where  n  is  the  number  of  grid  points. 

The  flow-chart  showing  the  solution  scheme  is  given  in  Figure  20. 
As  seen  in  the  flow-chart,  when  and  T2  reach  0.9  times  the 
governing  equations  are  modified. 

Program  Logic,  Description  of  the  Computing  Procedure 

For  each  set  of  given  physical  input  data  including  the  nominal 
film  thickness  h^ ,  the  program  calculates  the  pressure  distribution  in 
the  inlet  zone  including  the  inlet  pressure  p^  (at  =  a,  i.e. ,  h  -  h^) 
The  nominal  film  thickness  depends  slightly  on  the  inlet  pressure  at 
high  inlet  pressures  and  this  dependence  decreases  as  the  inlet  pressure 
increases.  Therefore,  h^  tends  to  an  asymptotic  value,  on  h^  versus  p^, 
curves  for  a  given  set  of  operating  conditions  (U,W,G).  The  first  step 
in  this  analysis  obtains  h^  versus  p^  curves  by  varying  h^  for  each  set 
of  physical  input  data.  Starting  h^  values  are  slightly  greater  than 
Grubin's  prediction  for  nominal  film  thickness,  and  decreased  by  small 
amounts,  until  the  conputer  program  "dumped"  because  of  the  logarithmic 
functions  in  Equations  (3-a)  and  (3-b).  It  is  clear  from  the  physical 
situation  that  decreasing  the  nominal  film  thickness  is  equivalent,  in 
a  sense,  to  increasing  pressure  (and  pressure  gradient). 
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Since  the  shear  stress  is  directly  related  to  the  pressure 
gradient,  these  increased  pressure  gradients  will,  eventually,  cause 
the  shear  stress  to  reach  the  limiting  value.  Therefore,  the  argument 
of  some  of  the  logarithmic  functions  in  Equations  (3-a)  and  (3-b) 
approach  zero.  The  ''dumping''  of  the  computer  program  dictates  a  limit 
to  the  decrease  of  h^,  vdiich  actually  comes  from  the  rheological  model 
as  and  ^2  approach  unity.  In  the  integrated  momentum  equation, 

Equation  (4A)  the  product  (p'h)  is  always  positive,  therefore  1X2!  is 
always  greater  than  |t^|.  A  physical  justification  for  this  is  to  run 
the  upper  surface  faster. 

Consequently,  it  is  expected  that  T2  will  approach  Xj^  before  x^^. 

In  order  to  handle  this  difficulty  in  the  numerical  solution  we  impose 
the  criteria  that  when  X2  ^0.9  Xj^  (and  x^  <  0.9  Xj^)  a  modified  version 
of  the  governing  equations  (3-a)  and  (3-b)  are  solved  where 
(1  -  t2)^  iln(l  -  X2)  ,  n  =  1,2  terms  are  omitted  based  on  the  limit  identity 


£im  X 
x^O 


n 


£n(x)  =  0 


If  then  x^  ^0.9  Xj^  (and  X2  ^  0.9  Xj^)  ,  we  impose  the  condition  =  0 

(i.e.,  p  constant),  as  required  by  the  integrated  momentum  equation, 

Equation  (4A) .  Imposing  p'  =  0  will  be  discussed  later  in  the  discussion 

section.  If  both  x^  and  X2  approach  Xj^,  then  p'  approaches  zero, 

provided  that  h  is  nonzero.  Although  the  numerical  problem  of  "dumping" 

due  to  the  logarithmic  functions  is  handled  in  this  manner  drastic  decreases 

in  h  still  do  not  occur, 
o 
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The  operating  value  o£  h  predicted  from  the  analysis  is  the  h 

°  dh 

where  >  100  MPa  (15  kpsi)  and  ^  approaches  zero  from  plots  like 
those  shown  in  Figures  2  and  3. 


30 


D.  Physical  Input  Data 

The  input  data  employed  in  the  analysis  are  shown  in  Tables 
I  through  IV.  The  lubricant  properties  are  representative  o£  typical 
lubricants  at  40C. 


Table  I.  Lubricant  and  Contact  Material  Properties 


2 

u  /mPas  (iVsBa/in  )  . 

410  (5.94  X  10"^) 

41  (5.94  X  10“^) 

a/GPa  ^  Cpsi 

32  (2.18  X  lO'"^) 

/I''lPa(psi) 

0 

6.9  (1000) 

0.69  (100) 

m/ dimensionless 

0.1 

0.05 

R/mm  (in.) 

12.7  (0.5) 

Contact  material 

Steel  on  steel 
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E.  Results 

h  versus  p -  curves 
0 _ 1 

Figure  21shows  the  effect  of  slide/roll  ratio,  2,  on  the  film 

thickness  inlet  pressure  relation,  for  the  physical  input  shown. 

For  large  h^  (i.e.,  small  p^) ,  change  of  I  does  not  have  much 

effect  on  the  inlet  pressure.  Although  it  is  not  shown  in  this  figure, 

the  solution  of  Reynold's  equation  with  a  Nev^donian  fluid  as  used  by 

Grub in  approaches  these  three  curves  at  low  inlet  pressure  values.  For 

inlet  pressures  around  130  MPa  (20  kpsi)  or  higher  this  relation  reduces 

to  h^  independent  of  p^.  A  straight  line  shown  in  Figure  21  represents 

this  Grub in  film  thickness.  At  different  inlet  pressure  ranges,  all 

three  curves  tend  to  go  to  an  asymptotic  h^  value.  As  shown  in  the 

figure,  increasing  E  reduces  the  film  thickness. 

In  Figure  22,  the  effect  of  limiting  shear  stress  properties  on 

film  thickness  is  shown.  The  physical  input  is  the  same  as  in  the 

previous  figure,  with  E  =  2  except  for  and  m  as  indicated.  The 

^o 

effect  of  increasing  and/or  m  is  qualitatively  the  same  as  the  effect 

o 

of  decreasing  E.  Xj^  and  m  can  be  varied  independently  but  their  effect 

o 

will  always  be  in  the  same  direction. 

The  Shear  Stress  Distribution 

The  shear  stress  distributions  in  the  inlet  zone  shown  in  Figure  23 
are  for  the  operating  conditions  represented  by  the  point  B  in  Figure  2. 
Notice  that  both  x^  and  X2  are  below  the  values  throughout  the  inlet 


zone. 
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In  the  case  o£  a  decreased  (i.e.,  increased  p)  the  shear 

stress  distribution  is  shown  in  Figure  24  for  the  conditions  repre- 

sented  by  point  A  in  Figure  21,  where  both  —  and  —  are  around  unity. 

This  point  is  chosen  as  the  point  defining  h^  predicted  by  our 

analysis  at  which  the  contact  will  operate  under  the  given  physical 

input.  The  pressure  gradient  (See  Figure  28)  is  zero  near  the  inlet 

(pressure  constant,  see  Figure  29) ,  and  h^  tends  to  go  to  an  asymptotic 

value  at  inlet  pressures  around  100  to  130  MPa  (15  to  20  kpsi) . 

The  shear  stress  distribution  under  pure  rolling  (point  C  of 

Figure  21)  is  shown  in  Figure  25.  Notice  the  signs  and  behavior  of 

and  T2  both  of  which  approach  zero  near  the  entrance  to  Hertz  contact 

resulting  in  a  zero  pressure  gradient  at  that  point. 

Pressure  Gradient  and  Pressure  Distribution 

Figure  26  represents  the  conditions  denoted  by  point  B  in  Figure  21 . 

The  related  pressure  distribution  is  given  in  the  Figure  27.  in  Figure  28, 

the  p'  distribution  for  the  conditons  represented  by  point  A  of  Figure  21 

is  sho™.  Notice  the  increase  in  the  p'  value  compared  with  that  in 

^max  ^ 

Figure  26  (although  not  by  order  of  magnitude  due  to  the  fact  that  h^ 
values  does  not  differ  much)  and  the  movement  of  the  peak  to  the  right. 
The  pressure  distribution  for  this  case  is  given  in  Figure  29 . 

Pressure  gradient  and  pressure  distributions  under  pure  rolling 
condition  are  shown  in  Figures  30  and  31,  for  the  conditions  represented 
by  the  point  C  of  Figure  21  (point  C,  where  p^  is  around  260  MPa,  is -out 
of  the  range  of  the  Figure  21) .  The  value  is  greater  than  the  one 

shown  in  Figure  26  by  one  order  of  magnitude .  The  related  pressure 
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distribution  is  shown  in  Figure  31.  This  resembles  Grub in's  infinite 

inlet  pressure  condition. 

H  versus  U,  and  H  versus  W  curves 
o  ’ _ o _ 

Figure  32  shows  the  dependency  of  on  U  for  w/L  =  87.6 

kN/m(500  Ib/in)  (p^  =  0.5  GPa  =  72  kpsi) ,  with  the  material  property 

G  =  6540,  and  for  two  different  slide/roll  ratios,  and  u  are  varied 

in  the  range  displayed  in  the  previous  tables  to  maintain  the  indicated 

U  values  [See  Table  IV].  The  limiting  shear  stress  parameters  are  kept 

constant.  Qualitatively,  since  the  effect  of  changing  the  limiting 

shear  stress  parameters  is  the  same  as  changing  the  slide -roll  ratio 

(See  Figures  21  and  2^  ,  these  variations  are  not  shown  on  versus  U 

curves.  The  conditions  for  Figure  33  differ  from  those  in  Figure  32 

only  in  that  the  load  is,  w/L  =  8.76  MN/m  (5  x  lo"^  Ib/in)  (or  Pj^  =  5 

GPa  =  727  kpsi) .  The  data  in  Figures  32  and  33  are  plotted  in  Figure  34 

to  show  the  dependence  of  the  film  thickness  on  load  which  is  only 

slightly  changed  from  the  classical  Grubin  result.  (Note  that  limiting 

shear  stress  parameters  =0.69  ^'lPa  m  =  0.05;  and  the  material 

o 

parameter  G  =  6540  remain  unchanged  in  Figures  32^  33  and  34.) 
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F.  Discussion 

The  basic  equations  employed:  Only  the  visco -plastic  portion 
o£  the  constitutive  equation,  Equation  2  was  used.  By  examining  the 
results ,  the  assumption  o£  omitting  the  elastic  term  o£  that  equation 
appears  to  be  justi£ied.  The  results  presented  permit  the  £ollowing 
order  or  magnitude  analysis: 

At  =  0.7  to  7  MPa  (10^  to  10^  psi) 

Ax^  -  0.25  to  0.025  mm  (lO'^  to  lO'^  in) 

u  =  s  0.25  to  2.5  m/s  (10  to  10  in/s) 

and  Tor  the  typical  pressure  range  encountered, 

G  =0.7  GPa  (10^  psi) 

oo 

T.  -  0.7  to  7  MPa  (10^  to  10^  psi) 

1j 

y  0.7  to  7  Pas  (10"^  to  10"^  Ibs/in^) 
as  —  approaches  unity  tanh  ^  \  will  be  at  least  in  the  range 

■^L  V^L/ 

o£  1  to  10  or  much  greater.  ThereTore,  i£  we  approximate  the  elastic 
term  in  the  constitutive  equation  as 

i_  dx  ^  _L  ^  i_  At  u 
G  Ht  G  dx  dt  G  AX^ 

OO  oo  1  00  1 
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0  3 

It  will  be  on  the  order  of  10  to  10  while  the  visco -plastic  term 
used  will  be 

or  more  as  approaches  unity. 

L 

Therefore  for  the  case  we  deal  with  in  this  analysis  the  elastic 
term  is  smaller  than  the  nonlinear  viscous  term  by  at  least  two  orders 
of  magnitude  due  to  the  fact  that  the  elastic  modulus  is  high,  and 
viscosity  is  low.  Consequently,  it  is  reasonable  to  omit  the  elastic 
term  of  the  constitutive  equation  in  this  inlet  analysis.  This 
assumption  may  not  be  valid  for  an  analysis  which  considers  the  Hertzian 
contact  zone  \diere  pressure  (therefore  viscosity)  is  much  higher  than 
it  is  at  the  inlet  zone. 

The  assunption  of  an  isothermal  flow  condition  can  be  justified 
for  the  rolling  case  by  using  the  results  of  Murch  and  Wilson  |2lJ  and 
Cheng  [22].  Their  thermal  reduction  factor  cj)^  is  available,  which  can 
be  multiplied  with  the  isothermal  calculation  of  h^,  to  find  the  actual  h^. 
For  our  high  viscosity  case  using  a  thermal  conductivity  of  mineral  oil 
0.12  W/m-C  (0.015  lb/°F-sec)  and  temperature  coefficient  of  0.04C"^ 

(0.02  F  ^) ,  (f)^  will  be  0.9  or  higher  where  u  =  200  in/sec  or  smaller. 

The  energy  dissipation  with  our  non-linear  fluid  will  be  less  than 
the  Nexrtonian  fluid  they  consider  for  other  conditions  equal.  Hence 
neglecting  the  inlet  zone  heating  effects  and  assuming  isothermal 
condition  is  valid  for  the  range  of  variables  used. 
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As  stated  in  the  above  discussion  o£  the  solution  technique, 
when  both  and  T2  approach  the  limiting  shear  stress  the  pressure 
gradient  was  assigned  a  value  of  zero.  As  seen  from  the  integrated 
momentum  equation  (Equation  4A)  ,  if  x^  and  X2  were  both  equal  to  the 
limiting  value,  the  pressure  gradient  would  have  to  be  zero  because  h 
is  nonzero.  However,  the  shear  constitutive  equation  (2-b)  shows  the 
shear  stresses  only  reach  the  limiting  value  as  the  shear  rate  approaches 
infinity.  Consequently  the  p'h  term  of  the  momentum  equation  must 
approach  zero  and  is  the  difference  between  two  large  nearly  equal  numbers. 
The  numerical  difficulties  this  causes  were  eliminated  by  sinply  forcing 
p'  to  zero  although  it  would  be  small  not  zero.  This  step  was  only 
required  very  near  the  Hertzian  inlet  zone  as  the  operating  film  thick¬ 
ness  was  approached,  (See  Figure  (28)).  when  p'  is  set  to  zero  the 
constitutive  equation  combined  with  momentum  Equation  (5A)  would  give 
a  constant  shear  rate  across  the  film,  hence  a  linear  velocity  profile, 
which  violates  continuity  because  h  is  not  constant  outside  the  Hertzian 
region.  The  numerical  convenience  of  forcing  p'  to  zero  under  these 
conditions  is  thought  to  have  no  major  effect  on  the  conclusions  of  the 
analysis  and  no  more  important  than  other  assumptions  inherent  in  a 

Grub in- type  solution  numerically  executed. 

The  effect  of  the  conventional  terms,  load,  speed  and  the  material 

parameters,  on  film  thickness  is  in  the  same  sense  as  expected  in  any 
Grub in  type  film  thickness  analysis.  For  a  fixed  material  parameter  G, 
the  dimensionless  film  thickness  H  depends  slightly  upon  the  load 
parameter  W  but  is  more  sensitive  to  the  changes  in  the  speed  parameter  U 
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as  seen  in  the  Figures  34,  33  and  32  respectively.  IVhen  the  sum 
velocity  approaches  zero  (or  the  vsicosity  y  approaching  zero) ,  the 
nominal  film  thickness  approaches  zero  (i.e. ,  dry  contact)  and  the  three 
lines  of  Figures  32  and  33  approach  each  other  for  smaller  film  thick¬ 
nesses  . 

The  effect  of  increasing  the  limiting  shear  stress  parameters , 

Tj^  and  m,  as  shown  in  Figure  32  is  to  cause  the  model  to  give  results 
o 

close  to  the  results  of  the  Newtonian  model.  When  either  or  m 

o 

(or  both)  is  increased  the  resulting  Xj^  value  at  that  particular  point 

is  increased,  hence  the  term  is  decreased  and  it  is  less  likely  to 

L 

approach  unity  with  increasing  |x|.  The  behavior  approaches  Newtonian 
viscous  behavior  and  the  asymptotic  tendency  of  h^  is  therefore  to  Grub in's 
predictions . 

The  Xj^  and  m  values  (at  a  specified  temperature)  are  available 
o 

from  the  experimental  our  previous  studies  [2].  Some  of  the  values 

used  in  this  analysis  are  approximated  from  the  data  on  the  lubricant 

5P4E  (polyphenyl  ether)  at  40C.  The  importance  of  the  slope  m  dominates 

the  importance  of  the  zero  pressure  value  because  near  the  inlet  to 

o 

the  Hertzian  contact  zone  the  pressure  increases  rapidly  and  x^  follows 
this  rapid  increase  through  the  parameter  m.  The  effect  of  x  is  far 

JL 

o 

out  in  the  inlet  zone  where  the  pressure  is  low  and  the  shear  stress  is 
low. 

X  . 

The  effect  of  increasing  sliding  speed  on  ,  i  =  1,2  is 

L 

qualitatively  the  same  as  the  effect  of  decreasing  the  limiting  shear 
stress  parameters.  Apparently  this  effect  is  not  altered  by  the  rolling 
speed  at  which  the  sliding  is  increased.  Decreasing  the  slide-roll  ratio 
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(or  increasing  the  limiting  shear  stress  parameters)  increases  the  h^ 
value  predicted  by  this  analysis,  and  shifts  this  asymptotic  value  to 
the  right  on  h^  versus  p^  curves.  This  is  due  to  the  fact  that  at  low 
slide/roll  ratios  and/or  with  high  limiting  shear  stress  parameters  the 
lubricant  behaves  like  a  Newtonian  fluid  and  therefore  at  the  asymptotic 
value  of  h^,  higher  pressures  (and  pressure  gradients)  are  possible  to 
achieve. 

The  effect  of  the  lubricant  viscosity  can  be  analyzed  in  the 
same  sense  as  the  effect  of  the  sliding  speed,  since  increase  in  either 
one  will  cause  the  shear  stress  to  increase.  It  was  observed  from  a 
few  solutions  of  the  model  that  for  materials  like  some  silicone  fluids 
with  extremely  high  values  and  low  values  the  shear  stresses  on 
the  surfaces  reach  the  limiting  value  even  at  unreasonably  high  h^ 
values  although  pressure  and  pressure  gradients  are  low.  With  this  type 
of  lubricant  our  assumption  of  omitting  the  elastic  term  from  the  con¬ 
stitutive  equation  may  not  be  valid.  The  high  viscosity  lubricants  work 
with  high  limiting  shear  stress  parameters  and  low  slide/roll  ratios, 
if  they  work  at  all. 
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G.  Conclusions 

This  analysis  confirms  the  idea  that  when  the  shear  stress 

reaches  the  limiting  value  in  the  inlet  zone  there  will  be  a 

reduction  in  the  nominal  film  thickness.  This  reduction,  for  most 

conventional  lubricants,  is  not  drastic.  For  the  slide/roll  ratio 

value  of  two,  it  is  found  that  the  nominal  film  thickness  is  about 

thirty  to  forty  percent  less  than  Grubin's  prediction  when  the 

-4 

dimensionless  film  thickness,  H,  is  in  the  order  of  10  ;  and  about 

twenty  percent  or  less  reduction  is  predicted  when  the  film  thickness 
is  in  the  order  of  10  ^  or  less . 

The  slide/roll  ratio,  and  the  lubricant  limiting  shear  stress 
are  the  two  newly  introduced  concepts  in  a  Grubin  type  film  thickness 
analysis.  Decreasing  the  slide/roll  ratio  and/or  increasing  the  limiting 
shear  stress  parameters  resulted  in  Newtonian -like  behavior.  This 
characteristic  behavior  is  implied  in  the  non-linear  visco-plastic 
rheological  model. 

Although  the  confirmation  of  the  existence  of  a  limiting  shear 
stress  is  relatively  new  in  the  field  of  tribology,  it  is  an  intrinsic 
material  property  of  the  lubricants  that  have  been  used  for  many  years. 
Therefore  it  is  reasonable  to  expect  as  concluded  that,  under  usual 
operating  conditions,  with  conventional  lubricants,  no  drastic  changes 
or  sudden  collapse  in  film  thickness  is  expected.  Hence  the  small 
(<  40  percent)  reduction  in  film  thickness  predicted  in  this  analysis 


is  reasonable. 
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APPENDIX  A 


Symbol : 


Source : 


Type; 


DESCRIPTION  OF  EXPERIMENTAL  FLUIDS 


Sun  Oil  Company 
Naphthenic  Base  Oil  R-620-15 


2 

Properties;  Viscosity  at  37. 8C,  m  /s 
Viscosity  at  98. 9C,  m^/s 
Viscosity  Index  (ASTM  D-2270) 
Flash  Point,  C 
Pour  Point,  C 
Density  at  20C,  Kg/m^ 

Average  Molecular  Weight 


24.1  X  10 


3.73  X  10' 


915.7 


Symbol ; 


Source : 


Type; 


MCS-460 


Monsanto  Company 
Synthetic  Hydrocarbon 


Properties;  Viscosity  at  37. 8C,  m  /s 

2 

Viscosity  at  98. 9C,  m  /s 
Viscosity  at  148. 9C,  m^/s 
Pour  Point,  C 
Density  25C,  Kg/m^ 


37.2  X  10 


4.0  X  10 
1.9  X  10" 


-29  to  -32 


932.7 
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Symbol : 

Santotrac  50 

Source : 

I^nsanto  Conpany 

Type: 

Synthetic  Cycloaliphatic  Hydrocarbon  Traction  Fluid 

Properties : 

2 

Viscosity  at  37. 8C,  m  /s 

34  X  10'^ 

Viscosity  at  98. 9C,  m^/s 

5.6  X  10  ^ 

Pour  Point,  C 

-37 

3 

Density  at  37. 8C,  Kg/m 

889 

Flash  Point,  C 

163 

Fire  Point,  C 

174 

Specific  Heat  at  37. 8C,  J/Kg'K 

2332 

Additive  package  includes:  Antiwear  (zinc  dialkyl 
dithiophosphate) ,  Oxidation  inhibitor.  Antifoam,  VI 

Improver  (Polymethacrylate) . 

Symbol : 

5P4E 

Type: 

Five -ring  Polyphenyl  Ether 

Source : 

Monsanto  Conpany 

Properties : 

2 

Viscosity  at  37. 8C,  m  /s 

363  X  10'^ 

2 

Viscosity  at  98. 9C,  m  /s 

13.1  X  10'^ 

Density  at  22. 2C,  Kg/m^ 

1205 

3 

Density  at  37. 8C,  Kg/m 

1190 

Flash  Point,  C 

288 

Pour  Point,  C 

4.4 
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Symbol : 
Type: 

Blend 

Properties : 


Source  of 
Polymer: 

Properties 
of  Polymer: 


Symbol : 
Type: 

Source : 
Properties : 


N3 

Blend  of  N1  and  Polyalkylmethacrylate 
(PL-4523) 

Viscosity  at  "hl.lC,  m^/s  182  x  10  ^ 

Viscosity  at  98. 8C,  m^/s  27  x  10  ^ 

Pressure  viscosity  coefficient 
(atmospheric  pressure  slope) 

at  37. 7C,  GPa'^  10.7 

at  98. 8C,  GPa"^  8.27 

Rohm  and  Haas  Company 

Viscosity  Average  Molecular  Weight  1.65  x  10^ 

Viscosity  m^/s  at  98. 9C  773  x  10  ^ 

Consists  of  19%  Polymer  in  solution 
with  a  paraffinic  hydrocarbon 


Krytox  143 -AB  (Lot  10) 
Perfluorinated  polyether 
DuPont  Con^any 

2 

Viscosity  at  37. 8C,  m  /s 
Viscosity  at  98, 9C,  m^/s 

3 

Density  at  24C,  kg/m 
Density  at  98. 9C,  kg/m^ 
V.I.  (ASTM  D-2270) 

Pour  point,  C 
Flammability 


96.6  X  10  " 

11.5  X  10'^ 

1890 

1760 

116 

-40 

does  not  bum 
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Symbol : 
Source : 
Type; 

Properties : 


XRM-177-F  (Lot  4) 

Mobil 

Synthetic  Paraffinic  Hydrocarbon 
Viscosity  at  37. 8C,  Pas 
Viscosity  at  98. 9C,  Pas 
Pour  Point,  C 
Density  37. 8C,  Kg/m^ 


376  X  10 
31.6  X  10 
<  -40 


838.9 
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where 


48 


Then  substituting  Equation  (4B)  into  Equation  (IB)  results  in 


^tanh-1 


i  dp 

dx. 


X3  +  F(x^) 


(SB) 


It  is  clear  from  Equation  (5B)  that  for  ^  =  0,  the  velocity  u^  is  a 

linear  function  of  x^,  which  with  the  no  slip  boundary  conditions 

^11  ^71 

gives  the  flowrate  Q  =  — f±  h^.  This  is  the  same  result  as  found 

with  the  usual  Newtonian  model  [14] .  Integrating  Equation  SB  from 
surface  1  with  respect  to  x^  gives  the  velocity  distribution 

u^Cx^)  - 


p'Xv  +  F'' 


£n 


p'Xv  +  F' 


p'x,  +  F' 


£n 


If  x^  =  h,  then  u^  =  u^j  and  it  will  be  the  Equation  (3-a)  of 
the  text. 

Integrating  the  velocity  distribution,  Equation  (6B) ,  in  con¬ 
tinuity,  Equation  (3B)  will  result  in  Equation  (3-b)  of  the  text. 


PRESSURE 


CONSTANT  PRESSURE 
SHEAR  STRESS 
APPARATUS 


Constant  Pressure  Stress -Strain  Apparatus,  1.1  GPa 
(Section  A:  modification  for  solid  samnle) . 
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\Sd'A  /  SS3dlS  dV3HS  ONIlllAin 


Limiting  Shear  Stress  vs .  Temperature  for  MCS  460 , 
Krytox,  Nl,  and  N1  +  4  percent  PAMA  (Cooling  History) 


LIMITING  SHEAR  STRESS  /  KPSI 


INSTRUMENT  HISTORY 


T  /  C 
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I _ I _ I - 

0.4  0  6  0.8 

PRESSURE  /GPa 


Figure  5.  Limiting  Shear  Stress  of  5P4E  (polyphenyl  ether)  vs, 
Pressure  for  two  Histories  and  Two  Temperature 


LIMITING  SHEAR  STRESS  /  KPSI 


history:  isothermal  COMPRESSION 

FROM  LIQUID 


Elastic  Shear  Ifodulus  o£  N1  (Naphthenic  Base  Oil)  vs 
Shear  Rate  at  Four  Pressures  and  21C 


SHEAR  STRESS  /  MF 


CELL  DEFLECTION  /  mm 


Figure  9. 


Shear  Stress -Cell  Deflection  for  Teflon  at  Four 
Pressures  and  21C 
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CELL  DEFLECTION  /  mm 

Figure  10.  Shear  Stress-Cell  Deflection  for  Acrylic  Plastic  at 
Three  Pressures  and  25C 


CELL  DEFLECTION  /  mm 


Figure  11.  Shear  Stress-Cell  Deflection  for  Nylon  at  Four 
Pressures  and  22C 


LIMITING  SHEAR  STRESS  /  KPSI 
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Figure  12 


•  PVC,  21  c 

■  TEFLON,  21  C 

★  Nl,  25C 

X  SANTOTRAC  50,  25  C 


i 
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0.5  1.0 


PRESSURE  /GPa 


Shear  Strength  or  Limiting  Stress  vs.  Pressure  for  Nl 
(Naphthenic  Base  Oil) ,  Santotrac  50  (Cycloaliphatic 
iiydrocarbon) ,  Polyvinyl  Chloride ,  and  Teflon 


Bdv\i  /  ssauis  av3HS 


Speed  of  0.22  in/s 
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Figure  14,  Traction  Curves  Predicted  from  Model  for  5P4E  at  45C,  Rolling  Speed 
of  0,22  m/s,  and  Hertz  Pressures  and  Side-Slip  Angle  as  Indicated, 


Compared  with  Measured  Data  from  Johnson  and  Tevaarwerk 
[81  for  indicated  Hertz  Pressure 
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Figure  16.  Traction  Coefficient  as  Predicted  from  Model  as  a  Function 
of  Slide-Roll  Ratio  and  Side-slip  Angle  in  radians  for  N1 
Mineral  Oil  at  30C,  1.0  m/s  Rolling  Speed,  0.86  CPa  Hertz 
Pressure  and  Spin  Radius  of  36  mm. 


Various  Relative  Side-slip  Angles,  N1  Mineral  Oil  at  30C, 
1.0  m/s  Rolling  Speed,  0.86  GPa  Hertz  Pressure  and  Constant 
Spin  Radius  of  36  mm 
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1.  Ti<0.9t^,  T2<0.9Tj^ 


2.  T-j  <  0. 9  T|^  ,  T2  >  0. 9  T|^ 


3.  T-,  >  0.9  T|  ,  ’''2  -  ^ 


Figure  20.  Solution  scheme 


Figure  21.  Film  thickaess  versus  inlet  pressure,  for  different  slide-roll 
ratios  fw/L  =  87.6  kN/m,  Pu  =  0.5  GPa,  u  =  2.54  m/s,  y  =  410 
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Figure  23.  Shear  stress  distribution  at  the  inlet  zone  for  the  conditions 
of  point  B  of  Figure  21 

★:  Limiting  shear  stress 

O :  Upper  surface  shear  stress 

• :  Lower  surface  shear  stress 


a 

Figure  24.  Shear  stress  distribution  at  the  inlet  zone  for  the  conditions 
of  point  A  of  Figure  21  (Symbols  as  in  Figure  23) 


SHEAR  STRESS  /MPa 
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Figure  25.  Shear  stress  distribution  at  the  inlet  zone  for  the  conditions 
of  point  C  of  Figure  21 (Symbols  as  in  Figure  23) 


PRESSU 


Figure  30.  Pressure  gradient  distribution  at  the  inlet  zone  for  the  conditions 
of  point  C  of  Figure  21 
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Figure  31.  Pressure  distribution  at  the  inlet  zone  for  the  conditions 
of  point  C  of  Figure  21 
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